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Abstract 

In this paper, Lie symmetry method is performed for the Hirota- 
Ramani (H-R) equation. We will find The symmetry group and optimal 
systems of Lie subalgebras. Furthermore, preliminary classification of its 
group invariant solutions, symmetry reduction and nonclassical symmeries 
are investigated. Finally the conservation laws of the H-R equation are 
presented. 
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1 Introduction 

In the present paper, we study the following equation 

U t -u x 2 t +au x (l-u t ) =0, (1) 

where a ^ and u(x, t) is the amplitude of relevant wave mode. This equation 
was introduced by Hirota and Ramani in [10]. Jie Ji obtained some travel- 
ling soliton solutions of this equation by using Exp- function method [12]. This 
equation is completely integrable by the inverse scattering method. Eq. (1) 
is studied in [10, 12, 9] where new kind of solutions were obtained. Hirota- 
Ramani equation is widely used in various branches of physics, such as plasma 
physics, fluid physics, quantum field theory. It also describes a variety of wave 
phenomena in plasma and solid state [10]. 

The theory of Lie symmetry groups of differential equations was developed by 
Sophus Lie [14], which was called classical Lie method. Nowadays, application 
of Lie transformations group theory for constructing the solutions of nonlinear 
partial differential equations (PDEs) can be regarded as one of the most active 
fields of research in the theory of nonlinear PDEs and applications. Such Lie 
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groups are invertible point transformations of both the dependent and inde- 
pendent variables of the differential equations. The symmetry group methods 
provide an ultimate arsenal for analysis of differential equations and is of great 
importance to understand and to construct solutions of differential equations. 
Several applications of Lie groups in the theory of differential equations were 
discussed in the literature, the most important ones are: reduction of order 
of ordinary differential equations, construction of invariant solutions, mapping 
solutions to other solutions and the detection of linearizing transformations for 
many other applications of Lie symmetries see [17, 5, 4]. 

The fact that symmetry reductions for many PDEs are unobtainable by 
applying the classical symmetry method, motivated the creation of several gen- 
eralizations of the classical Lie group method for symmetry reductions. The 
nonclassical symmetry method of reduction was devised originally by Bluman 
and Cole in 1969 [6], to find new exact solutions of the heat equation. The 
description of the method is presented in [7, 13]. Many authors have used the 
nonclassical method to solve PDEs. In [8] Clarkson and Mansfield have pro- 
posed an algorithm for calculating the determining equations associated to the 
nonclassical method. A new procedure for finding nonclassical symmetries has 
been proposed by Bila and Niesen in [1]. 

Many PDEs in the applied sciences and engineering are continuity equations 
which express conservation of mass, momentum, energy, or electric charge. Such 
equations occur in, e.g., fluid mechanics, particle and quantum physics, plasma 
physics, elasticity, gas dynamics, electromagnetism, magneto-hydro-dynamics, 
nonlinear optics, etc. In the study of PDEs, conservation laws are important 
for investigating integrability and linearization mappings and for establishing 
existence and uniqueness of solutions. They are also used in the analysis of 
stability and global behavior of solutions [2, 3, 19, 20]. 

This work is organized as follows. In section 2 we recall some results needed 
to construct Lie point symmetries of a given system of differential equations. 
In section 3, we give the general form of a infinitesimal generator admitted by 
Eq. (1) and find transformed solutions. In Section 4, we construct the optimal 
system of one-dimensional subalgebras. Lie invariants, similarity reduced equa- 
tions and differential invariants corresponding to the infinitesimal symmetries of 
Eq. (1) are obtained in section 5 and 6. Section 7, is devoted to the nonclassi- 
cal symmetries of the H-R model, symmetries generated when a supplementary 
condition, the invariance surface condition, is imposed. Finally in last section, 
the conservation laws of the Eq. (1) are obtained. 

2 Method of Lie Symmetries 

In this section, we recall the general procedure for determining symmetries for 
any system of partial differential equations see [17, 15, 5, 4]. To begin, let us 
consider the general case of a nonlinear system E of partial differential equations 
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of order n in p independent and q dependent variables is given as a system of 
equations 

A„(a;,uW)=0 ) v = !,-■■, I, (2) 

involving x — (x 1 , • • • , x p ) , u — (u 1 , ■ ■ ■ , u q ) and the derivatives of u with respect 
to x up to n, where represents all the derivatives of u of all orders from 
to n. We consider a one-parameter Lie group of infinitesimal transformations 
acting on the independent and dependent variables of the system (2) 

x l = x l + s£ l (x,u) + 0(s 2 ), i = l---,p, 

(3) 

v? = u 3 '+ sifP{x, u) + 0(s ), j = l---,q, 

where s is the parameter of the transformation and £ 8 , rf are the infinitesimals of 
the transformations for the independent and dependent variables, respectively. 
The infinitesimal generator v associated with the above group of transformations 
can be written as 

P Q 

v = j2e(x,u)d xi + j2<p a (x, U )d u «. (4) 

i— 1 a— 1 

A symmetry of a differential equation is a transformation which maps solutions 
of the equation to other solutions. The invariance of the system (2) under the 
infinitesimal transformations leads to the invariance conditions (Theorem 2.36 
of [17]) 

Pr (n) v[A„(xy n) )] = 0, i/ = l, whenever A 1/ (x,w (n) ) = 0, (5) 

where Pr'™^ is called the n th order prolongation of the infinitesimal generator 
given by 

Pr(")v = v + 22^(a;,«W)a t ,5, (6) 

Q = l J 

where J = ■ ■ ■ ,jk), 1 < jk < P, 1 < k < n and the sum is over all J's of 
order < #J < n. If # J = k, the coefficient ipj of d u ^ will only depend on 
fc-th and lower order derivatives of u, and 

v>i( X , «W) - r>j( Va - ^ ro + E (7) 

i=l t=l 

where uf := du a /dx l and u"^ := duj/dx 1 . 

One of the most important properties of these infinitesimal symmetries is 
that they form a Lie algebra under the usual Lie bracket. 
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3 Lie symmetries of the H-R equation 

We consider the one parameter Lie group of infinitesimal transformations on 

(a; 1 = x, x 2 = t 7 u = u), 

x = x + s£(x,t,u) + 0(s 2 ), 

t = x + sr](x,t,u) + C»(s 2 ), (8) 
u = x + stp(x, t, u) + 0(s 2 ), 

where s is the group parameter and £ = £, £ 2 = r\ and ip 1 = ip are the in- 
finitesimals of the transformations for the independent and dependent variables, 
respectively. The associated vector field is of the form: 

v = £(x, t, u)d x + r}(x, t, u)d t + <p(x, t, u)d u . (9) 

and, by (6) its third prolongation is 

Pr (3) v = v + <p x 8 Ux + tp* d Ut + if? d Ux2 + <p xt 8 Uxt + </ d Ut2 

+<(?' d Ux3 + V xH d u ^ t + <f* d Uxt2 + f 8 Ut3 . (10) 

where, for instance by (7) we have 



<p x = D x (ip-^u x -r]Ut)+^u x 2+r)u xt , 
Lp l = D t ((p-£u x -7iu t )+£u xt + riu t 2, 



V = D*(ip-£u x -riu t )+£u xt 3+r)u t 4. 



(11) 



where D x and D t are the total derivatives with respect to x and t respectively. 
By (5) the vector field v generates a one parameter symmetry group of Eq. (1) 
if and only if 

Pr (3) v[w t - u x 2 t + au x (l - u t )] = 0, 

(12) 

whenever u t — u x 2 t + au x (l — u t ) — 0. 
The condition (12) is equivalent to 

(1 - au x )Lp l + a(l - u t )ip x - f 1 = 0, 

(13) 

whenever u t — u x 2 t + au x (l — u t ) = 0. 

Substituting (11) into (13), and equating the coefficients of the various mono- 
mials in partial derivatives with respect to x and various power of u, we can 
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find the determining equations for the symmetry group of the Eq. (1). Solving 
this equations, we get the following forms of the coefficient functions 

—x 2t u 2x 

C=— C1 + C3, 77 = Cl t + c 2 , ip = (— + -- —)c 1 +c 4 . (14) 

where c 4 , c 2 , c 3 and c 4 arc arbitrary constant. Thus, the Lie algebra q of 
infinitesimal symmetry of the Eq. (1) is spanned bye the four vector fields 

1 2x 

vi=d x , v 2 = d t , v 3 = -d u , V4 = 3tdt xd x + (2t + u )d u . (15) 

a a 

The commutation relations between these vector fields are given in the Table 
1. The Lie algebra g is solvable, because if = (v,, [vj,Vj]) = [q,q], we have 
qW = ( Vl , • • • v 4 ), and fl ( 2) = [0 (1) , (1) ] = (-Vi - 2v 3 ,3v 2 + 2av 3 ,v 3 ), so, we 
have a chain of ideals D D 0. 



Table 1: The commutator table 





Vl 


v 2 


V3 


v 4 


Vl 











-vi — 2v 3 


v 2 











3v2 + 2av 3 


V3 











V3 


V 4 


vi + 2v 3 


— 3v2 — 2av 3 


-V3 






To obtain the group transformation which is generated by the infinitesimal 
generators v, for i = 1, 2, 3, 4 we need to solve the three systems of first order 
ordinary differential equations 



dx(s) 

ds 
dt(s) 

ds 
du(s) 

ds 



= £i(x(s),t(s),u(s)), x(0)=x, 
= r}i(x(s),i(s),u(s)), t(0)=t, 
= tpi(x(s),i(s),u(s)), u(0)=u. 



(16) 



Exponentiating the infinitesimal symmetries of Eq. (1), we get the one-parameter 
groups Gi(s) generated by Vj for i = 1, • • • , 4 



(17) 





(t,x 


u) - 


->■ (or + s,t, tt), 


G 2 


(t,x 


u) - 


-> (or,t + s,u), 


G 3 


(t,x 


u) - 


-> (x,t,u + s/a), 


G 4 


(t,x 


u) 1 


-> (xe~ s ,te 3s ,te 3s + 



1. 



Consequently, 
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Theorem 3.1 If u = f(x,t) is a solution of Eq. (1), so are the functions 
Gi (a) •/(*,*) = f(x-s,t), 

G 2 (s)-f(x,t) = f(x,t-s), (18) 
G 3 (s)-f(x,t) = f(x,t)+s/a, 

G 4 (s)-f(x,t) = e s f(xe s ,te- 3s ) + ^-(l-c 2s )+t(l-c- 2s ). 



4 Optimal system of the H-R equation 

In general, to each s-parameter subgroup H of the full symmetry group G of 
a system of differential equations in p > s independent variables, there will 
correspond a family of group-invariant solutions. Since there are almost always 
an infinite number of such subgroups, it is not usually feasible to list all possible 
group-invariant solutions to the system. We need an effective, systematic means 
of classifying these solutions, leading to an "optimal system" of group-invariant 
solutions from which every other such solution can be derived. [17] 

Definition 4.1 Let G be a Lie group with Lie algebra g. An optimal sys- 
tem of s— parameter subgroups is a list of conjugacy inequivalent s— parameter 
subalgcbras with the property that any other subgroup is conjugate to pre- 
cisely one subgroup in the list. Similarly, a list of s— parameter subalgebras 
forms an optimal system if every s— parameter subalgebra of g is equivalent to 
a unique member of the list under some element of the adjoint representation: 
l) = Ad( 5 (f))).[17] 

Theorem 4.2 Let H and H be connected s- dimensional Lie subgroups of the 
Lie group G with corresponding Lie subalgebras f) and f) of the Lie algebra g of 
G. Then H—gHg^ 1 are conjugate subgroups if and only if \) — Ad(g(f))) are 
conjugate subalgebras. [17] 

By theorem (4.2), the problem of finding an optimal system of subgroups 
is equivalent to that of finding an optimal system of subalgebras. For one- 
dimensional subalgebras, this classification problem is essentially the same as 
the problem of classifying the orbits of the adjoint representation, since each one- 
dimensional subalgebra is determined by nonzero vector in g. This problem is 
attacked by the naive approach of taking a general element V in g and subjecting 
it to various adjoint transformation so as to "simplify" it as much as possible. 
Thus we will deal with the construction of the optimal system of subalgebras of 
0- 
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To compute the adjoint representation, we use the Lie series 

e 2 

Ad(exp(evi)vj) = v, - e[v h vj] + y [v;, [v i; vj]} , (19) 

where [vj,Vj] is the commutator for the Lie algebra, e is a parameter, and 
i, j = l, 2, 3, 4. Then we have the Table 2. 



Tabic 2: Adjoint representation table of the infinitesimal generators 



Ad 


Vl 


v 2 


V3 


v 4 


Vl 


Vl 


v 2 


V3 


v 4 + e(vi + 2v 3 ) 


v 2 


Vl 


v 2 


v 3 


v 4 — e(3v 2 + 2av 3 ) 


V3 


Vl 


v 2 


V3 


v 4 — ev 3 


V 4 


vi - e(vi + 2v 3 ) 


v 2 + e(3v 2 + 2av 3 ) 


v 3 + ev 3 


v 4 



Theorem 4.3 An optimal system of one- dimensional Lie algebras of the H-R 
equation is provided by 

1) v 4 , 2) avi+^v 2 + v 3 , 3)avi+v 2 , 4) Vi 



Proof: Consider the symmetry algebra g of the equation (1) whose adjoint 
representation was determined in table 2 and let F? : g — V g defined by v n> 
Ad(exp(evi)v) is a linear map, for i = 1, • • • ,4. The matrices Mf of Ff, i = 
1, • • • , 4, with respect to basis {vi, • • • , v 4 } are 



Ml = 



1 

1 












— e 


-2s 
1 



,M| = 



1 














1 





3e 








1 


las 











1 



M\ = 



c 


— e 




1 









1 

1 £ 

1 



-3 s 
-2e 



(e--l) 







Let V = 5Z i=1 ajVi is a nonzero vector field in g. We will simplify as many 
of the coefficients at as possible by acting these matrices on a vector field V 
alternatively. 

Suppose first that a 4 ^ 0, scaling V if necessary we can assume that a 4 = 1, 
then we can make the coefficients of Vi, v 2 and V3 vanish by Mf and Mf . And 
V reduced to case 1. 

If 04 = and a 3 ^ 0, then we can not make vanish the coefficients of vi 
and v 2 by acting any matrices Mf. Scaling V if necessary, we can assume that 
a 3 = 1 and V reduced to case 2. 
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If <i4 = <23 = and a2 ^ 0, then we can not make vanish the coefficient of 
vi. Scaling V if necessary, we can assume that ai = 1 and V reduced to case 3. 

The remaining one-dimensional subalgebras are spanned by vectors of the 
above form with = 03 = 02 = 0. If a\ ^ 0, we scale to make a\ = 1, and V 
reduced to case 3. □ 



5 Symmetry reduction of the H-R equation 

Lie-group method is applicable to both linear and non-linear partial differential 
equations, which leads to similarity variables that may be used to reduce the 
number of independent variables in partial differential equations. By determin- 
ing the transformation group under which a given partial differential equation 
is invariant, we can obtain information about the invariants and symmetries of 
that equation. 

Symmetry group method will be applied to the (1) to be connected directly to 
some order differential equations. To do this, a particular linear combinations of 
infinitesimals are considered and their corresponding invariants are determined. 
The equation (1) is expressed in the coordinates (x,t,u), so to reduce this 
equation is to search for its form in specific coordinates. Those coordinates 
will be constructed by searching for independent invariants (y, v) corresponding 
to the infinitesimal generator. So using the chain rule, the expression of the 
equation in the new coordinate allows us to the reduced equation. Here we will 
obtain some invariant solutions with respect to symmetries. First we obtain the 
similarity variables for each term of the Lie algebra 0, then we use this method 
to reduced the PDE and find the invariant solutions. 

We can now compute the invariants associated with the symmetry operators, 
they can be obtained by integrating the characteristic equations. For example 
for the operator v 2 + v 3 = d t + \d u characteristic equation is 

dx dt adu 

T ~ T ~ ~~ T' ( ' 

The corresponding invariants are y = x, v — u — ^ therefore, a solution of our 
equation in this case is u = v(y) + The derivatives of u are given in terms of 
y and v as 

U X =Vy, U x 2=Vyy, 14^=0, U t = ~ ■ (21) 

Substituting (21) into the Eq. (1), we obtain the ordinary differential equation 
(a — l)v y + l/a = 0, the solution of this equation is v — ^ryy+c. Consequently, 
we obtain that 

T t 

+ -+c. (22) 



a(l — a) a 
All results are coming in the tables 3 and 4 
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Table 3: Reduction of Eq. (1) 



operator 


y 




V 


u 


Vl 


t 




u 


v{y) 


v 2 


X 






v(y) 


V4 


xt 1 / 3 


(u — 2x/a 


)t -l/ 3+t 2/3 


v(y)t 1/3 + 2x/a - t 


vi + av3 


t 


u 


— X 


v(y) +x 


v 2 + v 3 


x 


u - 


-t/a 


v(y) + t/a 


Vl + v 2 


X — t 




u 


v(y)/a 



Table 4: Reduced equations corresponding to infinitesimal symmetries 



operator similarity reduced equations 

Vi v y = 

V2 av y = 

V 4 t- 2/3 V + (Xt- 1/3 - 3)Vy + Xt- 2/3 Vyy + (tt^^Vy + 2) (6 - ' 3 V - t" 1 / 3 ^) = 3 

vi + av3 (1 — a)Vy + a — 

V2+V3 (a — l)v y + l/a = 

Vl + V2 — Vy + Vyyy + aVy(l + Vy ) = 



6 Characterization of differential invariants 

Differential invariants help us to find general systems of differential equations 
which admit a prescribed symmetry group. One say, if G is a symmetry group 
for a system of PDEs with functionally differential invariants, then, the system 
can be rewritten in terms of differential invariants. For finding the differential 
invariants of the Eq. (1) up to order 2, we should solve the following systems of 
PDEs: 

dl_ dl_ ldl_ 3t d£_ x dl (2t u _^_ ) ^l m) 

dx' dt' adu dt dx a du 1 

where / is a smooth function of (x, t, u) 1 

dh dh Idh 3t dh _ x 8h + ... + (2 _ 2ut) ^L (24) 
dx ' dt a du' dt dx du t ' 

where I\ is a smooth function of (x, t, u, u x , u t ) 1 

dl 2 dl 2 1 dl 2 „ dl 2 dl 2 dl 2 r dl 2 

-57, 3t—-x—^ u xt - 5ti«- — , (25) 

dx dt a du dt dx du xt du u 

where h is a smooth function of (x,t 7 u,u x ,u t ,u XX7 u x t 7 utt)- The solutions of 
PDEs systems (23), (24) and (25) coming in table 5, where * and ** are refer to 
ordinary invariants and first order differential invariants respectively. 
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Table 5: differential invariants 



vector field 


ordinary invariant 


1st order diff. invariant 


2nd order diff. invariant 


Vl 
V2 
V3 
V 4 


t, u 
X, u 
X, t 

tx 3 , - t + u)x 


*, U x ,Ut 
*, U x ,Ut 
*, U x ,Ut 

* r 2 v - i "'TT 1 

*i x "a: a ' ^ 


*,**, U xx ,U xt ,U t t 
*, **, U xx ,U x t, Utt 

*,**,u xx ,u xt ,u u 
* ** r s ii ^i- H*t 



7 Nonclassical symmetries of the H-R equation 

In this section we would like to apply the nonclassical method to the H-R equa- 
tion. The graph of a solution 

u a = f a (x u ---,x n ), a = l,---,q (26) 

to the system (2) defines an p-dimensional submanifold r / C R p x R q of the 
space of independent and dependent variables. The solution will be invariant 
under the one-parameter subgroup generated by vector (4) if and only if is 
an invariant submanifold of this group. By applying the well known criterion 
of invariance of a submanifold under a vector field we get that (26) is invariant 
under vector (4) if and only if / satisfies the first order system Eq of partial 
differential equations 

p 

Q a (x,u,uW) = ^ a \x,u)-J2e(x,u)u? -0, a=l, •••,<? (27) 

i=l 

known as the invariant surface conditions. The q-tuplc Q = (Q 1 ,-''!*? 9 ) is 
known as the characteristic of the vector field (4). In what follows, the n-th 
prolongation of the invariant surface conditions (27) will be denoted by Eq\ 
which is a n-th order system of partial differential equations obtained by ap- 
pending to (27) its partial derivatives with respect to the independent variables 
of orders j < n — 1. 

For the system (2), (27) to be compatible, the n-th prolongation Pr^v of 

(n) 

the vector field v must be tangent to the intersection E n Eq 

PrWv(A v )| £nB c ; , ) =0, v = l,.-.,l. (28) 

If the equations (28) are satisfied, then the vector field (28) is called a nonclas- 
sical infinitesimal symmetry of the system (2). The relations (28) are gener- 
alizations of the relations (5) for the vector fields of the infinitesimal classical 
symmetries. A similar procedure is applicable to the case of the nonclassical in- 
finitesimal symmetries with an evident difference that in general one has fewer 
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determining equations than in the classical case. Therefore, we expect that 
nonclassical symmetries are much more numerous than classical ones, since any 
classical symmetry is clearly a nonclassical one. The important feature of de- 
termining equations for nonclassical symmetries is that they are nonlinear, this 
implies that the space of nonclassical symmetries does not, in general, form a 
vector space. For more theoretical background see [18, 1]. 
Consider the system E of second order equations 

ut - v t + au x (l - u t ) = 0, u xx - v = (29) 

obtained from the H-R equation. If we assume that the coefficient of d t of the 
vector field (4) does not identically equal zero, then for the vector field 

v = £(x, t, u, v)d x + d t + ip(x, t, u, v)d u + ip(x, t, u, v)d v (30) 

the invariant surface conditions are 



u t + £u x = ip, v t + £v x = ip (31) 

The equations (28) take the forms 

(2 (u t - l/2)u 2 x a - u t (u x - v x ))£ u - u t i\) u - i/j t 
+ ((1 - 2 u t )au x + u t )f u + (u x (l - u t )a - Ut)ip v 

- l)a 2 + 2 (u t - l/2)u 2 x a - u t (u x - v x ))£ v (32) 
+ (u 2 x (u t - l)a 2 + ((1 - 2 u t )u x - v x (u t - l))a + u t )f v 
+ (1 - au x )(f t - a(u t - l)<f x + u x (u t - l)a£, x + {u 2 x a - u x + v x )£ t = 0, 

and 

^) / u> x £ tXX 2 u x ^ xu -\- u x f uu / u> x £ tUU 4" u xx f u 2 u xx ^ x -f- v x tf vv 
J r v xx fv ~\~ 2 u x f xu + 2 v x ip xv — 2 u x v x £ xv + 2 u x v x ip uv 2 u x v x £ uv (33) 
3u xx u x £u ^u xx v x ^ v u x v x ^ vv v xx u x ^ v + f xx 0. 

After inserting ip and its derivatives, as determined by the equation (33), in 
to (32) and substituting v x — u xxx ,v xx = u xxxx , and equating the coefficients 
of the various monomials in partial derivatives with respect to x and various 
power of u, we can find the determining equations. Solving this equations, we 
get four Algebraic equations equal to zero. This means that no supplementary 
symmetries, of non-classical type, are specific for our model. 

Now assume that the coefficient of dt in (30) equals zero and try to find the 
infinitesimal nonclassical symmetries of the form 

v = d x + if(x, t, u, v)d u + ip(x, t, u, v)d v (34) 

for which the invariant surface conditions are the following ones 

u x = if, v x =ip (35) 
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Relations (28) lead to the system of equations for the functions ip and tp 

{u 2 x {u t - l)a 2 + {{-v x - 2u x )u t +v x + u x )a + u t )tp v 

+ {{u x - 2 u x u t )a + u t )ip u + (u x (u t - l)a - u t )tp v (36) 

+ (1 - au x )ip t - a(u t - l)ip x -ip t - Utipu = 0, 

and 

—i> + Vxx + 2 u x tp xu + 2 tp xv v x + u 2 x ip uu + 2 u x v x ip uv . . 

+u xx tp u + V 2 x lf vv + V xx (fi v =0. 

Similar the previous case, we can find determining equations. Solving this equa- 
tions, we get the following form of the coefficient functions 

ip = c\u V t(l - 2u t )c\ + c 2 (38) 

where ci and ci are arbitrary constant. So the system (29) admits the classical 
symmetry Vi and nonclassical symmetry v 2 : 

X 

vi = d x + d u , v 2 = d x - ('■ u-t + 2 tu t )d u (39) 

a 

8 Conservation laws of the H-R equation 

Many methods for dealing with the conservation laws are derived, such as the 
method based on the Noether's theorem, the multiplier method, by the rela- 
tionship between the conserved vector of a PDE and the Lie-Bcklund symmetry 
generators of the PDE, the direct method, etc. [17, 2, 3, 19]. 

Now, we derive the conservation laws from the multiplier method. 



Definition 8.1 A local conservation law of the PDE system (2) is a divergence 
expression 

Di&[u] - D^iu] + ■■■ + D n $ n [u] = (40) 

holding for all solutions of the system (2). In (40), <fr l [w] = $ l (x, it, d u , ■ ■ ■ , <9£), 
i = 1, • • • ,n, are called fluxes of the conservation law, and the highest-order 
derivative (r) present in the fluxes is called the order of a conservation 

law. [3] 

In particular, a set of multipliers {A„[[/]}(, =1 = {A v (x, U, du, • • • , dlj)} v=1 
yields a divergence expression for the system A v (x, u^) such that if the identity 

A V [U]A V [U] = Di&[U] (41) 

holds identically for arbitrary functions U(x). Then on the solutions U(x) = 
u(x) of the system (2), if A„[C/] is non-singular, one has local conservation law 
A„[u]A„[u] = Di&lu] =0. 
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Definition 8.2 The Euler operator with respect to U J is the operator defined 

by 

^=^7-A^- + - + (-irA 1 -A s ^ + - (42) 



li- 



fer j = 1, • • • ,q. [3] 



Theorem 8.3 The equations E u jF(x, U, du, ■ ■ ■ , dfj) = 0, j = 1, • • • , q hold 
for arbitrary U (x) if and only if F(x, U, du, ■ • • , <9j>) = D i ^! t (x, U, du, ■•■ , d^ 1 ) 
holds for some functions ^ % {x, U,du,- ■ ■ , <9^ _1 ), i = 1, • • • <?• [3] □ 

Theorem 8.4 A set of non-singular local multipliers {A I/ (x, U,djj, - ■ ■ , ^[})}L=i 
yields a local conservation law for the system A^(x, u( n )) if and only if the set 
of identities 

E w (K v {x, U,d ur --, d^A^x, u^)) =0,j = l,...q, (43) 

holds for arbitrary functions U(x). [3] □ 

The set of equations (43) yields the set of linear determining equations to find 
all sets of local conservation law multipliers of the system (2). Now, we seek all 
local conservation law multipliers of the form A = £(x, t, u, u x , u t , u xx , u xt , u tt ) 
of the equation (1). The determining equations (43) become 

EuMx,t, U, U x , U t , U xx , U xt , U tt )(U t - XJ xH + aU x (l - U t ))] = 0, (44) 

where U (x, t) are arbitrary function. Equation (44) split with respect to third 
order derivatives of U to yield the determining PDE system whose solutions are 
the sets of local multipliers of all nontrivial local conservation laws of second 
order of H-R equation. 

The solution of the determining system (44) given by 

ci U xx + \c 2 {2 tU tt + U t -l) + c 3 u tt7 (45) 

where Ci, C2 and C3 are arbitrary constant. So local multipliers given by 

l)S = U xx , 2)£ = U tt , 3) ^ = tU tt + ±U t - \, (46) 

Each of the local multipliers £ determines a nontrivial two-order local conserva- 
tion law D t *S> + D x § = with the characteristic form 

A* + Dx* = m - U xH + aU x (l - U t )), (47) 

To calculate the conserved quantities ^ and $, we need to invert the total 
divergence operator. This requires the integration (by parts) of an expression 
in multi-dimensions involving arbitrary functions and its derivatives, which is 
a difficult and cumbersome task. The homotopy operator [20] is a powerful 
algorithmic tool (explicit formula) that originates from homological algebra and 
variational bi-complexes. 
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Definition 8.5 The 2-dimcnsional homotopy operator is a vector operator 
with two components, (#^ t) /, ^u( x ,t)f) > wnere 

<U/=f(E^/)[A«]f -nd <U/-f(Ee ) /)[A.]^.(48) 
The x-integrand, /, is given by 

Mf M| kl -l k 2 

fc 1= lfc 2 =0 ii=0i 2 =0 ° U x k it k 2 

where M\ , M| are the order of / in u to x and t respectively and combinatorial 
coefficient 

(h + «2\ Ai + - h - ii — 1\ 
£?( x ) = B(ii,i 2 , ki,k 2 ) = {n ){ / „ ^r: 1 - 1 J . (50) 



Ai + 



fcl J 



Similarly, the t-integrand, /*• /, is defined as 



Mj" hi fe 2 -l 

$/= E E (E E^w^k-^- 1 )^-, (5i) 

fe 1= 0fc 2 = l 11=0 12=0 OU x k lt k 2 

where B^fo, ii, fe, k\). 

For instance we apply homotopy operator to find conserved quantities >J> and 
$ which yield of multiplier £ = u tt . We have 

/=u M (u t -u x 2 t + au x (l-u t )), (52) 

the integrands (49) and (51) are 

fijf= auu t i - auu t u t 2 - \uu xi? + \u t u xt 2 + \u x u t3 - \u xt u t 2, 



ft]f = v.* - u t u x 2 t + au x u t - au x u\ + \uu x 2 t 2 - auu xt + auu xt u t 

+ \u x U xt 2 - ^U x 2U t 2, 

apply (48) to the integrands (53), therefore 

^u(x,t)f= \^uu t 2 - \auu t u t 2 - \uu xt z + \u t u xt 2 + \u x u t z - \u xt u t 2, 

^u\x,t)f = h u t - \ u t^ x 2 t + \au x u t - \au x u1 + \uu x 2 t 2 - \auu xt 
+\auu xt u t + \u x u xt 2 - \u x 2U t 2, 



(53) 



(54) 
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so, we have the first conservation low of the H-R equation respect to multiplier 



D x (\auu t 2 - \auu t u t 2 - \uu xt s + ^u t u xt 2 + \u x u t z - \u xt u t 2) 

+D t (\v% - \u t u x 2 t + \au x u t - \au x u\ + \uu x 2 t 2 - \auu xt 
+ \auu xt u t + \u x u xt 2 - \u x 2u t 2) = 0. 

Similarly, conservation law respect to multiplier £ = u xx is 

D x (\u x u t + \au x - \au x u t - huu xt + \auu x u xt ) 
+D t (\uu xx - \auu x u t 2 - \u xx ) = 0. 

and conservation law respect to multiplier £ = tu u + \u t — \ is 

Ar( — \ au + \auu t + ^atuu t 2 — ^atuu t u t 2 — \auv% — \uu xt 2 
-Uuu xt s + ±tu t u xt 2 - -^u t u xt + \u x u t 2 + \tu x u t 3 + \u xt 
-\tu xt u t 2) = 0. 



+D t ( - \u + \uu x 2 t + ±tUU x 2 t 2 + \tU x U xt 2 + J^U x U xt + i 

-j^u x 2u t - \atuu xt + \atuu t u xt + \tu\ - \tu t u x 2 t 
— \tu x 2u t 2 + ^atu x ut — \atu x uf) = 0. 



(55) 



(56) 



(57) 
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